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^ , In this talk we discuss the process e+e — — > it near threshold. In particular we 

discuss a quark mass definition which is a generalization of the static PS mass. 

The new definition allows us to calculate recoil corrections to the static PS mass. 
C _ ~) ' Using this result we calculate the cross section of e + e — — > tt near threshold at 

Q\ ' NNLO accuracy adopting three alternative approaches, namely (1) fixing the pole 

CD , mass, (2) fixing the PS mass, and (3) fixing the new mass which we call the PS 

f^*^ ■ mass. We demonstrate that perturbativc predictions for the cross section become 

f"^ ■ much more stable if we use the PS or the PS mass for the calculations. A careful 

analysis suggests that the top quark mass can be extracted from a threshold scan 

at NLC with an accuracy of about 100 - 200 MeV. 
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1 Introduction 

1.1 On NLC, FMC, LEG and Tevatron 

The top quark physics will be one of the main subject of future e + e~ and 
/j + /i~ colliders such as the Next Linear Collider (NLC) and the Future Muon 
Collider (FMC). The goals are to measure and to determine the properties 
of the top quark which was first discovered at the Tevatron El with a mass of 
m = 174.3 ± 5 GeVa. Although the top quark will be studied at the LHC and 
the Tevatron (RUN-II) with an expected accuracy for the mass of 2 — 3 GeV, 
the most accurate measurement of the mass with an accuracy of 0.1% (100 — 
200 MeV) is expected to be obtained only at the NLC0. 
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1.2 Why do we need 10 accuracy in the top mass? 

The top quark mass enters the relation between the electroweak precision ob- 
servables indirectly through loops effects. The global electro-weak fit of the 
Standard Model requires to have very accurate input data in order to make 
a constraint for the masses of undiscovered particles, such as the Higgs boson 
or other particles. The increase in the accuracy of the top quark mass will 
improve the limits on the Higgs mass. In addition we can study possible devi- 
ations from the Standard Model through anomalous couplings, CP violation, 
or extra dimensions. 

1.3 LO and NLO cross section 

Due to the large top quark width, the top-antitop pair cannot hadronize into 
toponium resonances. The cross section appears therefore to have a smooth 
line-shape showing only a modfirate IS peak. In addition the top quark wid 
serves as an infrared cutoffura and as a natural smearing over the energy 
As a result _thc nonperturbative QCD effects induced by the gluon condensate 
are small US, allowing us to calculate the cross section with high accuracy by 
using perturbative QCD even in the threshold region. Many theoretical styj 
at LO and NLO have been donp-ip-the past for the total cross section at 
for the-jBiOBjicntum distribution 119'tll, also accounting for electro- weak correc- 
tions liaH'EfLjajpjd-fjor the complete NLO correction including non-factorizable 
correctionsll3'll3 ! t3'lla. It was proven tr^pthojomufactorizable corrections cancel 
in the inclusive cross section at NLOy'BBli3BB. 

1.4 On NNLO results 



Recently, the NNL Q-^n^l^ 
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nclusive threshold production cross section 
The results of the NNLO analysis are 
summarized in a review article K3. To summarize the results for a standard 
approach using the pole mass, the NNLO corrections are uncomfortably large, 
spoiling the possibility for the top quark mass extraction at NLC with good 
accuracy because the IS peak is shifted by about 0.5 GeV by the NNLO, the 
last known correction. One of the main reasons for this is the usage of the pole 
mass in the calculations. It was realized that such type of instability is cs 
by the fact that the pole mass is a badly defined object within full QCDt 

In this talk we discuss a definition of the quark mass alternative to the 
pole_mass. It is the so-called potential subtracted (PS) mass, suggested in 
Ref.l23. In contrast to the pole mass the PS mass is not sensitive to the non- 
perturbative QCD effects. We derive recoil corrections to the relation of the 



pole mass to the PS mass and demonstrate that perturbative predictions for 
the cross section become much more stable at higher orders of QCD (shifts are 
below 0.1 GeV) if we use the PS mass for the calculations. This understanding 
removes one of the obstacles for the accurate top quark mass measurement 
and it can be expected that the top quark mass will be extracted from a 
threshold scan at NLC with an accuracy of about 100 — 200 MeV. We have 
to mention that the necessity of isolating the IR contributions in the mass 
calculatio n a n d t he co nsideration of a short distance mass have been studied 
intensivelyEJoOaCJ. The applica^n&MJthiLPS, LS and IS* mass ham been 
reported recently by several groups lH3'E£rEZH23'EZI and reviewed in Ref. lj. The 
detailed comparison of our results using the PS mass with results of other 
groups have been performed in Ref. E3 so that we refer the reader to this 
reference for details. 

In this talk we review Ref. E3. In particular, we discuss a quark mass 
definition (we call it PS pass) which is a generalization of the static PS mass 
proposed by M. Beneketia. The new definition allows us to calculate recoil 
corrections to the static PS mass. Second, we calculate the cross section for 
e + e~ — ► ti near threshold with NNLO accuracy using three alternative mass 
schemes, i.e. (1) the pole mass, (2) the static PS mass, and (3) the new ES 
mass. The results of the first approach has already been reported in Ref. c3, 
the results for the top quark pauL.pxoduc.tion near threshold using the static 
PS mass were presented in Refs.cjOEaEa, the third alternative schemes were 
discussed in Refs.Ea. 



2 Top quark production near threshold 
with NNLO accuracy 

In this section we consider the cross section of the process e + e~ — > tt in the 
near threshold region where the velocity v of the top quark is small. It is 
well known that the conventional perturbative expansion does not work in the 
non-relativistic region because of the presence of the Coulomb singularities 
at small velocities v — > 0. The terms proportional to (a s /v) n appear due to 
the instantaneous Coulomb interaction between the top and the antitop quark. 
The standard technique for re-summing the terms (a s /v) n consists in using the 
Schmdinger equation for the Coulomb potential and in finding the Green func- 
tion uB. The Green function is then related to the cross section by the optical 
theorem. In order to determine NLO corrections to the cross section we need 
to know the short-distance correction to the vector currentLJ, the NLO correc- 
tion to the .Coulomb potential c2l, and the contribution of the non-factorizable 
corrections 0Jll3'LZl. It was proven that the non-factorizable corrections cancel 



in the inclusive cross section at NLO and beyond 0000-00. At NNLO 
the situation is more complicated. One of the obstacles for a straightforward 
calculation are the UV divergences coming from relativistic corrections to the 
Coulomb potential. This problem can be solved by a proper factorization of the 
amplitudes and by employing effective theories. We want to sketch the deriva- 
tion of the inclusive cross section. The inclusive cross section can be obtained 
from the correlation function of two vector currents jy.(x) = t(x) r y IJ .t(x), 

IVb 2 ) = iJd 4 xe^(0\T{j^x), ]u (0)}\0). (1) 

The first step is to expand the Lagrangian and the currents 

ji = h^ = a^aiX - —^^a l (iT>) 2 x (2) 

om z 

consistently in 1/m. The useful language for treating the one and two noa.- 
relativistic quark(s) is provided by the NRQCD ffl and the PNRQCD E3, 
respectively. After the expansion, the cross section reads 



R = <j{e + e — > tt)/a p t 
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where a pt — Ana 2 /3s, eQ is the electric charge of the top quark, N c is the 
number of colors, -y/i = 2m + E is the total center-of-mass energy of the quark- 
antiquark system, m is the top quark pole mass and T is the top quark width. 
The unknown coefficient C{tq) can be fixed by using a direct QCD calculation 
of the vector vertex at NNLO in the so-called intermediate region eScS and by 
using the direct matching procedure suggested in Ref.E3. 

The function G(f, f'\E + iT) is the non-relativistic Green function, ft 
satisfies the Schrodinger equation 

(H - E- iT)G(r, f'\E + iT) = S(r- f '), (4) 

H is the -OQiLjEalativistic Hamiltonian of the heavy quark-antiquark system, ft 
is shownl^SEjEJ that the Schrodinger equation can be reduced to the equation 

(H 1 -E 1 )G 1 (r,r'\E 1 )=5 3 (r-r') (5) 

with the energy E\ — E + E 2 /4m, E = E + iT, and with the Hamiltonian 

Hi = h Vc(r) + — Vo(r) 
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r 
Vc(r) = V (r) I 1 + as ^l° F (2A) ln(A) + ai) (6) 
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/3 2 (^41n 2 ( M V) + y ) + 2(ft + 2/3bo x ) ln(/*'r) + a 2 j 



where // = /ie 7E , /i is the renormalization scale, and 75 is Euler's constant. 
The QCD beta function coefficients /3o and /?i and the coefficients ai and a 2 
are listed later this talk. The final expression for the NNLO cross section is 
given by 
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with Gi(r ,ro|£^i) being the solution of Eq. (13) at_c. = r' = r . For the 
numerical solution we use the program derived in Rcf.E.3 by one of the authors. 

3 On the mass definitions 

The top quark mass is an input parameter of the Standard Model. Although 
it is widely accepted that the quark masses are generated due to the Higgs 
mechanism, the value of the mass cannot be calculated from the Standard 
Model. Instead, quark masses have to be determined from the comparison of 
theoretical predictions and experimental data. 

It is important to stress that there is no unique definition of the quark mass. 
Because the quark cannot be observed as a free particle like the electron, the 
quark mass is a purely theoretical notion and depends on the concept adopted 
for its definition. The best known definitions are the pole mass and the MS 
mass. However, both definitions are not adequate for the analysis of top quark 
production near threshold. The pole mass should not be used because it has 
the renormalon ambiguity and cannot . Jac determined more accurately than 
300 - 400MeVBB (see also Refs. 230-13). Actually, we may relate the pole 
mass with some short distance mass like the MS mass M which for the "large 
/?o approximation" (see e.g. Ref.Ej) is given by 

m pole = Mh + Y,m n (^) ]. (8) 



The Borel image of the pole mass to po 1c reads 

TOp le(u) = Y^ m n — r (9) 



The function m po i e (M) ~ T(l — 2u) has poles at real and positive values u 
1/2, 3/2, . . .. Therefore, the inverse Borel transformation 



™ P oic = / duexp ( j m po ie(w) (10) 

generates ambiguities of the order 

<5m«exp(- JwAqcd- (11) 

\Poa s J 

These ambiguities_api known in the literature as renormalon ambiguities (see 
for example Refs.EJ'tiil). 

The MS mass is an Euclidean mass, defined at high virtuality, and therefore 
destroys the non-relativistic expansion. Instead, it_was recently suggested to 
use threshold masses like the low scale (LS) massEil, the potential subtracted 
(PS) mass E3, or one half of the perturbative mass of a fictious l 3 5*i ground 
state_icalled IS* mass) El In this talk we focus on the PS mass suggested in 
Ref.B, 



i r w< i*t d 3 k 

2 J W) 

where Vc is the quark-antiquark Coulomb potential 



mps = m polc - <5mp S with Sm PS = -- / ^-^ Vb(|fc|) (12) 



3.1 On renormalons and soft QCD effects cancellation 

In order to understand why this mass is better defined than the pole mass and 
to see that the pole mass is very sensitive to long distance effects, it is enough 
to consider the one-loop expression for the self energy diagram. Taking the 
residue in ko, one obtains a soft self energy contribution which comes from 
momenta k with \k\ < Aqcd, 

6m = Ava s C F J J^^ 2 = ^A QCD . (13) 

|*:|<Aqcd 

We observe that the pole mass has a non-perturbative uncertainty of order 
Aqcd which then penetrates into consequent perturbative QCD calculations. 



This uncertainty or ambiguity is caused by the poles of the Borcl transforma- 
tion as shown before. However, it is easy to realize that the PS mass is free of 
this ambiguity. Indeed, the term 8m in Eq. (|13|) and therefore the renormalon 
ambiguities associated with it cancel in the definition of the PS mass as given 
in Eq. ([L2|) as well as in the combination 2m po i e + V(r) which is known as 
static energy. We can see this by applying the Borel technique also to the 
QCD potential 

d 3 k ( AnC F a s {k 2 )\ ,,. 



The Borel image of the potential V(r) = ^ V n (f5Qa s /ATi) n in the "large j3q 
limit" is given by 






u 



V{u)=^V n — . (15) 

with the function V(u) ~ T(l — 2u) having again poles at real and positive 
values u — 1/2, 3/2, . . .. The inverse Borel transformation 

V(r) = fduexp (4^-1 V(u). (16) 

J \/3 a s J 

generates ambiguities of the order 

5V(r) « exp ( -^- ) w Aqcd- (17) 

The coefficient in V(u) differs from m po i e (w) by a factor —2. Therefore we ob- 
serve t he can cellation of the renormalon ambiguity in the combination 2m po i e + 
F(r)B@>@. 

The definition in Eq. JO) has been given in Ref. c3 but has akeady been 
discussed implicitly in Ref IS. The remarkable step made in Ref.Ea is to use 
this definition beyond one-loop order. It has been proven in Ref. 123 that the 



cancellation of the infrared QCD contributions to the PS mass in Eq. (12) 
holds even at higher loop orders. Still the definition in Eq. ( |l2| ) is valid only in 
the static approximation and does not contain \/m corrections like the other 
threshold mass definitions. We therefore suggest to extend the definition given 
in Ref. El. 

3.2 The PS mass definition 

Our objective is that the definition should be gauge independent and well- 
defined within quantum field theory so that radiative and relativistic correc- 
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Figure 1: The general structure of the self energy diagram of a quark 

tions can be calculated in a systematic way. Our definition is given by 

"ips = m po i e - (5mpg with Sm-^ = £ soft {$) , (18) 

where £ S oft is the soft part of the heavy quark self energy which is defined as 
the part where at least one of the heavy quark propagators is on-shell. Starting 
point for our considerations is therefore the self energy of an on-shell quark 
with mass m and momentum p (i.e. p 2 = to 2 ) as shown in Fig. |l| and written 
down as 

M 






n 



{-ig sl a "T an ) (19) 



where the last factor is a non-commutative product with decreasing index n. 
The line momenta k n are linear combinations of the gluon loop momenta l m , 
the particular representation is specified by the structure S. The symbol {l m } 
means the set of all these loop momenta, the same symbol is used for the 
Lorentz and color indices. In general we have M < N which means that line 
momenta can be correlated. The momenta of the virtual quark states are given 
byp n =p + k n . 

The quark is considered to be at rest, p — (to, 0), and it interacts with a 
number of gluons. The subdiagram S displayed in Fig. fil describes the inter- 
action between the gluons. In general the quark lines between the interaction 
points represent virtual quark states. However, if the virtual quark comes very 
close to the mass shell and the total momentum of the cloud of virtual gluons 



becomes soft, this situation gives rise to long-distance nonperturbative QCD 
interactions. The described (virtual) contributions results in the soft part of 
the self energy, E so ft . So we define 



N „ M 
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One can derive this expression from Eq. (y_9|) by using the identity 

a V ■ = -^b 2 - m 2 ) + P f -^-^) (21) 

yr — to z + ze \p z — m z / 

and the fact that the principal value integral does not give any infrared sensitive 
contribution. The delta function can be used to remove the integration over 
the zero component of ki . In order to parameterize the softness of the gluon 
cloud we impose a cutoff on the spatial component, \ki\ < /!/, and indicate 
this by a label /i/ written at the upper limit of the three-dimensional integral. 
This cutoff fj,f is also known as factorization scale. So we can rewrite Eq. (20) 



^B0ft(^ fif) = -^J2j J2^ V ^ P) 



(22) 
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.M-l , 4 
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The range of the index m is reduced by one which indicates that one of the 
loop momenta is extracted as line momentum of the i-th line. In the following 
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Figure 2: Leading order contribution to the quark self energy (a) and to the quark-antiquark 
potential (b). The cross indicates the point where we cut the quark line by imposing an on- 
shell condition to the virtual quark state. The gluon propagator can be decomposed in a 
Coulomb propagator (c) and a transverse propagator (d). 



we deal with the different realizations of this compact expression. As we will 
see explicitly, the function V(k,p) occurring as integrand can be seen as quark- 
antiquark potential where we have summed over the spin of the tensor product 
of a final state with an initial state. .Because the static quark-antiquark po- 
tential is used in a similar way in Ref.tEa, we recover the result of Ref.tZl in the 
static limit. 

3.3 The one-loop contribution 

The leading order contribution to the self energy of the quark is given by 

=»> = * / ^-">^ T ->jrhi { - i!hrT ' , i (24) 

where Feynman gauge is used for the gluon. The soft part of it is given by 

1 f»f d 3 k 
2~J (2^) 

where 



Esoft(AM/)= ■■— I 77^^y{k,p), V(k,p)=V+(k,p)+V_(k,p) (25) 



v ± (k, P ) = — g ^ (Vw2 + K2T2w L (26) 

2mym 2 + K?{ym 2 + n 2 — m) 

and k = \k\. The procedure of taking the soft part by setting a cut is shown 
in Fig. |2J(a-b). The two potential parts are known as the scattering potential 
V + (k,p) and the annihilation potential V_(fc,p) and correspond to the two 
zeros k± — ±\/n 2 + m 2 — m of the delta function. Integrating the potential 
up to the factorization scale (if, we obtain 



Z7r l m V rn I m V m z ' 
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Figure 3: Two-loop contributions to the quark self energy 

The expansion of this expression in small values of Hf/rn results in 

The first term reproduces the result given in Ref. Ej to leading order in a s 
while the second term is the recoil correction to the static limit in this order of 
perturbation theory. This second term is related to the Breit-Fermi potential 
but does not coincide with it. 



3-4 Two loop contributions 

To take a step beyond the leading order perturbation theory, we consider two- 
loop diagrams for the heavy quark self energy as shown in Fig. Q We calculate 
them in Coulomb gauge, even though we stress that our final result is gauge 
invariant. The gluon propagator in Coulomb gauge is given by 

jXab -rah / k-h\ 

G&(k) = ^, G<${k) = ^ (^ - ^ J , ij = 1,2,3. (29) 

The use of Coulomb gauge splits up the gluon propagators into a Coulomb term 
(Coulomb gluon) and a transverse term (transverse gluon) where the first one 
couples to the quark via the time components only. This splitting is shown in 
Fig.|(b-d). 

In cutting the quark line in all possible ways we obtain a lot of diagrams 
from the ones shown in Fig. 0. However, we find that the final contribution of 
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the two abelian diagrams in Fig. 0(a-b) to the soft part of the self energy is 
suppressed by //?/m 2 . Therefore it turns out that the only abelian diagrams 
which can give a non-suppressed contribution to the soft part of the quark self 
energy are the diagrams containing the vacuum polarization of the gluon as 
shown in Fig. [|(d-f). The simple calculation of these diagrams within the MS 
scheme, accounting only for light fermion loops, gluon loop (and ghost loop if 
Feynman gauge is used) results after renormalization in 
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1 f d 3 k I 4:TrC F a s (fi) 
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(30) 

This result has been anticipated because the expression in the curly brackets 
of the integrand reproduces the next-to-leading order correction to the QCD 
Coulomb potential. 

For the only non-abelian diagram in Fig. Bhc) the use of Coulomb gauge 
leads to seven diagrams. However, direct calculations show that only the dia- 
gram where the transverse gluon is joined by Coulomb gluons on both sides of 
the quark line contributes to the order g^/if/m. For this diagram and therefore 
for the whole non-abelian contribution up to this order we obtain 

Esoft - 8m »*' (31) 

This result has been anticipated, too, to be minus one half of the non-abelian 
correction to the QCELCoulomb potential, which is known in the literature 
(see for example Refs.l£3'E5), 

,na_ 1 r d 3 k J Tr 2 a 2 s C F C A \ _ o? s C F C A 2 



^~— 2 ] (2^ j m\k\ J ~ ^^ fif - (32) 

3. 5 Our final result 

Summarizing all contribution up to NNLO accuracy, we obtain 
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where m is the pole mass, fi is the renormalization scale, /if is the factorization 
scale, and 

Co = 1, C = 0, C = — — , 



d = 01 - 2/3 ln f^ , Cj = C A y, 

C 2 = a 2 - 2(2ai/3 + /?i) (in ( ^ 



+4/? 2 (ln 2 (^)-21n(^)+2). (34) 



The constants a\ and a 2 are given by 

31 r 20 t m 
ai = ~^Ca - —T F N F , 



EBB 



/4343 
fl2= lT62- 



, 7T 4 22 \ , /1798 56 \ 
4^ 2 - — + -C 3 ) C\ - ^ — + — C 3 J C A T F N F 



20„ „ \ 2 /55 



Q T F N F I - ( y - 16Ca I C F T F N F , (35) 

the coefficients of the beta function are known as 

0o = y Ca - I n f t f> 0i = yCl " y CaTfNf ~ ^ c fTfN f (36) 

where C F — 4/3, Ca — 3 and T F — 1/2 are color factors and N F = 5 is 
the number of light flavors. The coefficients C\ and C 2 in Eq. ( J33| ) have been 
derived in Ref. E3 by using known corrections to the QCD potential. In this 
work we have derived the coefficients Cg, C' ' , and C[. Note that our result 
can be represented in a condensed form as 

1 f- lf d 3 k / -. - \ 

m ws Qi,)-m = --J j^{Vc(\k\) + V R (\k\) + V NA (\k\)) (37) 

where the first term Vc is the static Coulomb potential, Vr is the relativistic 
correction (which is related to Breit-Fermi potential but does not coincide 
with it), and Vna is the non-abelian correction. Note that we did not include 
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electroweak corrections neither in this mass relation nor in the cross section 
we will look on later. Finally (which is actually not our result) the relation 
between the pole mass and the MS mass is given by the three-loop relationL3E3 

mpol ° -1 + | (—) + (— Y (-1.0414JV> + 13.4434) 



m(m) 3 V 7r / V 7T 

/ QJ \ 3 

+ [—) (0.65277VJ, - 26. 655N F + 190.595) (38) 

where a s — a s (jn) is taken at the MS mass. 

Using the relations between the masses, we fix the MS mass to take the 
values m(m) — 160 GeV, 165 GeV, and 170 GeV and use Eq. ( |38| ) to determine 
the pole mass at LO, NLO, and NNLO. This pole mass is then used as input 
parameter m in Eq. ( p3[ ) to determine the PS and PS masses at LO, NLO, 
and NNLO. The obtained values for the PS and PS mass differ only in NNLO. 
The results of these calculations .axa-collected in Table flL together with the 
estimates for "large /3o" accuracy ejEj. Note thai, the same values for the MS 
mass have been used for Tables 2 and 3 in Ref. Ej. The obtained mass values 
can now be used for the analysis in the following section. 

Taking the PS mass instead of the PS mass, we observe an improvement 
of the convergence. The differences for the mass values in going from LO to 
NLO to NNLO to the "large 0o" estimate for e.g. m(m) — 165 GeV read 
7.64 GeV, 1.64 GeV, 0.52 GeV, and 0.25 GeV for the pole mass, 6.72 GeV, 
1.21 GeV, 0.29 GeV, and 0.08 GeV for the_PS mass and finally 6.72 GeV, 
1.21 GeV, 0.27 GeV, and 0.08 GeV for the PS mass. 

In Fig. |j we show the difference between the PS and the PS mass (in 
GeV) as a function of the factorization scale [if (solid line) at /i = 15 GeV. It 
is interesting to observe that the non-abelian part of the difference between the 
PS and the PS mass (dotted line in Fig. |j) can be as large as 200 MeV. But 
the recoil correction cancel in part the non-abelian one and therefore the final 
difference is not more than 50 MeV. The dependence of m-pa — tops on the 
renormalization scale at [if — 30 GeV is given in Fig. |] by the dashed curve. 

4 NNLO results 

4-1 The scheme with the pole mass 

The top quark cross section at LO, NLO, and NNLO is shown in Fig. || as 
a function of the center-of-mass energy. For the top quark pole mass we use 
m t = 175.05 GeV, for the top quark-.width T t = 1.43 GeV, and for the QCD 
coupling constant a s (mz) — 0.119EI3. Different values [i = 15 GeV, 30 GeV, 
and 60 GeV for the renormalization scale are selected because they roughly 
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167.97 
167.95 


168.05 
168.03 


167.44 


169.05 


169.56 


169.80 


165.0 


171.72 


172.93 


173.22 
173.20 


173.30 
173.28 


172.64 


174.28 


174.80 


175.05 


170.0 


176.92 


178.17 


178.47 
178.45 


178.55 
178.53 


177.84 


179.52 


180.05 


180.30 



Table 1: Top quark mass relations for the MS, PS, PS, and the pole mass at LO, NLO, 
NNLO, and "large /3q" accuracy in GeV. We fix the MS mass to be m(m) = 160 GeV, 
165 GeV, and 170GeV and find the pole mass at LO, NLO, and NNLO from the three-lorm 
relation in Eq. (j38j). The PS and PS masses are derived from the pole mass by using Eq. (B3j) 
(without or with the 1/m contributions, respectively). We use the QCD coupling constant 
Ois(mz) = 0.119, (J, = m(m), and fif = 20 GeV for the factorization scale in the PS and PS 



correspond to the typical spatial momenta for the top quark. For solving the 
Schrodinger equation we use the program written by one of the authors lj. 
Note that we do not take into account an initial photon radiation which would 
change the shape of the cross section. This can be easily included in the Monte 
Carlo simulation. 

The NNLO curve modifies the line shape by the amount of 20 — 30% which 
is as large as the NLO correction. It also shifts the position of the IS peak by 
approximately 600 MeV. These large shifts of the peak positipn were expected. 
As we discussed above (and is well-known in the literature E3Eil), the pole mass 
definition suffers from the renormalon ambiguity. The top quark pole mass 
cannot be defined better than O(Aqcd)- Large NNLO corrections and a large 
shift of the IS resonance can spoil the top quark mass measurement at the 
NLO 



4-2 The scheme with the PS mass 

In this subsection we discuss the calculation scheme using the PS mass. We 
redefine the pole mass through the PS mass using the relation given in Eq. (33) 
without the 1/m contributions and then use the PS mass as an input parameter 
for our numerical analysis at LO, NLO, and NNLO. In Fig. g we show the top 
quark cross section expressed in terms of mps(/i/) at LO, NLO, and NNLO 
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(like in Fig. pi) as a function of the center-of-mass energy. We take mps (/i/ = 
20GeV) = 173.30 GeV which corresponds to the pole mass m — 175.05 GeV 
according to Table [y. In looking at Fig. g we observe an improvement in the 
stability of the position of the first peak in comparison to the previous analysis 
as we go from LO to NLO to NNLO. Actually, for the three values fj, = 15 GeV, 
30 GeV, and 60 GeV we obtain the maxima of the IS peak for NNLO at s max = 
347.32 GeV, 347.41 GeV, and 347.48 GeV while the maximal values are given 
by -Rmax = 1.379, 1.184, and 1.088, respectively. This demonstrated that a 
large variation in the renormalization scale [i gives rise only to a shift of about 
160 MeV for the 15 peak position at NNLO while the variation for i? max is 
still large. 

4-3 The scheme with the PS mass 

In this subsection we discuss the calculation scheme where we use the PS mass. 
We redefine the pole mass through the PS mass by using the relation given in 
Eq. (p3|) and then use the PS mass as an input parameter for the numerical 
analysis at LO, NLO, and NNLO. In Fig. [?] we show the top quark cross 
section expressed in terms of nipg(^j) at LO, NLO, and NNLO (like in Fig. ||) 
as a function of the center-of-mass energy. We take Wpg(/xj = 20 GeV) = 
173.28 GeV. Again we observe a very good stability of the position of the first 
peak as we go from LO to NLO to NNLO, similar to the one observed for the 
PS mass case. Studying the size of the NNLO corrections to the peak positions 
we conclude that the current theoretical uncertainty of the determination of 
the PS mass from the IS peak position is about 100 — 200 MeV. 

5 Conclusion and discussions 

We discussed the so-called potential subtracted (PS) mass suggested in Ref.Ea 
as a definition of the quark mass alternative to the pole mass. In contrast to 
the pole mass, this mass is not sensitive to the non-perturbative QCD effects. 
We have derived recoil corrections to the relation of the pole mass with the 
PS mass. The main result for this is Eq. (|33|). In addition, we demonstrated 
that, if we use the PS or the PS mass in the calculations, the perturbative 
predictions for the cross section become much more stable at higher orders of 
QCD (shifts are below 0.1 GeV). This understanding removes one of the ob- 
stacles for an accurate top mass measurement and one can expect that the top 
quark mass will be extracted from a threshold scan at NLC with an accuracy 
of about 100 -200 MeV. 
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Figure 4: The difference between the PS and the PS mass (in GeV) as a function of the 
factorization scale fif (solid line) at /j, = 15 GeV. The dotted line shows only the non-abclian 
part of the difference. The dependence of «}pg — m PS as a function of the normalization 
scale /i at /if = 30 GeV is shown as dashed line. 
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Figure 5: The scheme using the pole mass: shown is the relative cross section R(e + e~ — > tt) 
as a function of the center-of-mass energy in GeV for the LO (dashed-dotted lines), NLO 
(dashed lines), and NNLO (solid lines) approximation. We take the value mt = 175.05 GeV 
for the pole mass of the top quark, T t = 1.43 GeV for the top quark width, a s (rriz) = 0.119 
and different values fj, = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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Figure 6: The scheme using the PS mass: shown is the relative cross section R(e + e~ — » tt) 
as a function of the center-of-mass energy in GeV for the LO (dashed-dotted lines), NLO 
(dashed lines), and NNLO (solid lines) approximation. We take the value mpg = 173.30 GeV 
for the PS mass of the top quark, Y t = 1.43 GeV for the top quark width, a s (mz) = 0.119 
and different values fj, = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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Figure 7: The scheme using the PS mass: shown is the relative cross section i?(e + e~ — > tt) 
as a function of the center-of-mass energy in GeV for the LO (dashed-dotted lines), NLO 
(dashed lines), and NNLO (solid lines) approximation. We take the value m-p^ = 173.28 GeV 

for the PS mass of the top quark, Ft = 1.43 GeV for the top quark width, u s (mz) = 0.119 
and different values fi = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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